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Abstract 

This paper considers the application of multiple-input multiple-output (MIMO) techniques to non-orthogonal 
multiple access (NOMA) systems. A new design of precoding and detection matrices for MIMO-NOMA is proposed 
and its performance is analyzed for the case with a fixed set of power allocation coefficients. To further improve 
the performance gap between MIMO-NOMA and conventional orthogonal multiple access schemes, user pairing 
is applied to NOMA and its impact on the system performance is characterized. More sophisticated choices of 
power allocation coefficients are also proposed to meet various quality of service requirements. Finally computer 
simulation results are provided to facilitate the performance evaluation of MIMO-NOMA and also demonstrate the 
accuracy of the developed analytical results. 

I. Introduction 

Non-orthogonal multiple access (NOMA) has recently received considerable attention as a promising 
enabling technique in fifth generation (5G) mobile networks because of its superior spectral efficiency Jl] 
and GO. The key idea of NOMA is to explore the power domain, which has not been used for multiple 
access (MA) in the previous generations of mobile networks. Specifically NOMA users in one cell are 
served by a base station (BS) at the same time/code/frequency channel, and their signals are multiplexed 
by using different power allocation coefficients. The novelty of NOMA comes from the fact that users 
with poorer channel conditions are allocated more transmission power. In this way, these users are able 
to decode their own messages by treating the others’ information as noise, since the power level of their 
messages is higher. On the other hand, the users with better channel conditions will use the successive 
interference cancellation (SIC) strategy, i.e., they first decode the messages to the users with poorer channel 
conditions and then decode their own by removing the other user’s information. 
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The concept of NOMA can be linked to many well-known methods used in previous communication 
systems. For example, NOMA downlink transmissions resemble Cover and Thomas’s description of 
broadcast channels provided in [j3|. Another example is that the use of SIC has been extensively 
investigated in conventional multiple-input multiple-output networks, particularly in V-BLAST systems 
0. The superimposed messages transmitted in NOMA systems also resemble the concept of hierarchical 
modulation widely used for digital video broadcasting IQ. But unlike these existing techniques, NOMA 
seeks to strike a balance between throughput and fairness. For example, the transmission power allocated 
to the users in NOMA systems is inversely proportional to their channel conditions, which is important to 
ensure that all the users are served simultaneously. On the other hand, conventional opportunistic schemes 
prefer to give more power to users with better channel conditions, which can improve the overall system 
throughput but deteriorate fairness. 

The impact of path loss on the performance of NOMA has been characterized in [[6]| by assuming that 
users are randomly deployed in a cell, which has demonstrated that NOMA can outperform conventional 
orthogonal multiple access (OMA) schemes. In Q the implementation of NOMA has been considered 
in a scenario with two base stations, and the design of uplink NOMA has been proposed in 0. The 
user fairness of NOMA has been considered in by studying the impact of different choices of power 
allocation coefficients. In iTlOl a cognitive radio inspired NOMA scheme was proposed, in which the power 
allocation coefficients are chosen to meet the predefined users’ quality of service (QoS) requirements. 

In this paper, we focus on the application of multiple-input multiple-output (MIMO) to NOMA 
downlink communication systems. The concept of MIMO-NOMA has been validated by using systematic 
implementation in lUTI and llT2l . which demonstrates that the use of MIMO can outperform conventional 
MIMO-OMA. Compared to these existing works, the contributions of this paper are as follows: 

• We first consider a general NOMA downlink scenario, in which all users participate in NOMA with 
a fixed set of power allocation coefficients. A new design of precoding and detection matrices is 
proposed, and the impact of this design on the performance of NOMA is characterized by using the 
criteria of outage probabilities and diversity orders. The provided analytical and numerical results 
demonstrate that MIMO-NOMA can achieve better outage performance than conventional MIMO- 
OMA, even for users that suffer strong co-channel interference. 

• To enlarge the performance gap between MIMO-NOMA and MIMO-OMA, user pairing is applied to 
NOMA. Analytical results, such as an exact expression for the average sum-rate gap between MIMO- 
NOMA and MIMO-OMA and its high SNR approximation, are developed. These analytical results 
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demonstrate that the design of user pairing in NOMA systems is very different from conventional user 
scheduling scenarios. Conventionally it is preferable to schedule the users whose channel conditions 
are superior, but in the context of NOMA, it is important to schedule users whose channel conditions 
are very distinct. This is consistent with the findings obtained for single-antenna NOMA cases in JT| 
and ifTOll . 

• Inspired by the concept of cognitive radio networks, more sophisticated choices for the power 
allocation coefficients are proposed. Particularly we consider two types of constraints for the power 
allocation coefficients. One is to meet a predefined QoS requirement, i.e., a user’s rate supported by 
NOMA is larger than a targeted data rate. The other is to meet a more dynamic QoS requirement, 
where the user’s rate supported by NOMA needs to be larger than that supported by conventional 
OMA. Analytical results are developed for both scenarios to facilitate performance evaluation. 


II. System Model With Fixed Power Allocation 

Consider a downlink communication scenario, where a BS equipped with M antennas communicates 
with multiple users equipped with N antennas each. To make the NOMA principle applicable to this 
scenario, the users are randomly grouped into M clusters with K users in each cluster. The signals 
transmitted by the BS are given by 


x = Ps, 


( 1 ) 


where the Mxl vector s is given by 


0Y,1S1,1 + ' 

• + ®l,KS\,K 
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Sm 


( 2 ) 


where s rn ± denotes the information bearing signal to be transmitted to the A’-th user in the m-th cluster, 
ctij denotes the NOMA power allocation coefficient, and the design of the M x M precoding matrix P 
will be discussed in the next section. 

Without loss of generality, we focus on the users in the first cluster. The observation at the k -th user 
in the first cluster is given by 


yr.fc = H 1)fc Ps + n 1)fe , (3) 

where H ljfe is the N x M Rayleigh fading channel matrix from the BS to the A-th user in the first cluster, 
and ni fc is an additive Gaussian noise vector. Denote by vi^ the detection vector used this user. After 
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applying this detection vector, the signal model can be rewritten as follows: 

v i,kYi,k = v^ fc Hi ife Ps + vf fc ni )fe . 

Denote the z-th column of P by p,. The above signal model can be rewritten as follows: 

v ? fc yi,fc = v u H i ,fcPl (a!l,lSl,l + • • • + ai,KSl,K) 

M 


(4) 

(5) 


+ V U- H lT-PmSm + vf fc ni ife . 

m =2 

The channel conditions are crucial to the implementation of NOMA. Without loss of generality, we assume 
that the effective channel gains are ordered as follows: 

Iv^Hi^pil 2 > > |v^Hi )X pi| 2 , (6) 

and following the principle of NOMA, the users’ power allocation coefficients are ordered as follows: 


« 1,1 < • • • < <Xi,k- 


In this section, constant power allocation coefficients will be considered, and more sophisticated choices 
will be used in Section [VJ It is worth pointing out that optimizing power allocation according to 
instantaneous channel conditions can be used to further improve the performance of MIMO-NOMA, 
which is beyond the scope of this paper. 

Based on the above signal model, the signal-to-interference-plus-noise (SINR) for the AT-th ordered 
user in the first cluster is given by 


SINR 1K = 


(7) 


l v Mf H l,ifPl| 2 «l. 


K 


|2 1 


E*=7 Kk^kpiI^I + E"=2 Iv^h^pJ 2 + K* 

where p denotes the transmit signal to noise ratio (SNR). 

The k -th user, 1 < k < K, needs to decode the messages to the users with poorer channel conditions 
first, before detecting its own. The messages S\j, K > j > (k + 1), will be detected at the /c-th user with 
the following SINR: 


SINR{ k = 


( 8 ) 


l v ffc H i,fcPi| 2 «ij 


121 


Ei=i l v ^ H rfcPil 2a L + E m=2 |vf fc H 1 , fcPm | 2 + |v 1)fc 

If the message s i j can be decoded successfully, i.e., log(l + SINR [ k ) > Rij, then it will be removed 
from the k -th user’s observation, where R. t j denotes the j-th user’s targeted data rate. This SIC will be 
carried out until the k -th user’s own message is decoded with the SINR, SINR\ k . 
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The first user in the first cluster needs to decode all the other users’ messages with SINR{ 1; K > j > 2. 
If successful, it will decode its own message with the following SINR: 


SINR \! = 


I I „ 12 2 

v l,l-Hl,lPl| a i, 


E !=2 l v ffl H l,lPm| 2 + l v l,! 


121 
1 P 


(9) 


The design of the precoding and detection matrices will be discussed in the following section. 


III. Design of Precoding and Detection Matrices 


To completely remove inter-cluster interference, the precoding and detection matrices need to satisfy 
the following constraints: 

v".H ./ P,„ = 0, (10) 

for any m ^ i. 

In order to reduce system overhead caused by acquiring channel state information (CSI) at the BS, it 


is assumed that the BS does not have the global CS 


which leads to the following choice of P: 


P = I 


Mi 


where I M is the M x M identity matrix. The above choice means that the BS broadcasts the users’ 
messages without manipulating them. The advantage of this choice is that it avoids asking the users to 
feedback all their CSI to the BS, which consumes significant system overhead. 

With this choice of P, the constraints on the detection matrices in (fTOl) become 



m,ik 


o, 


(ii) 


where h m ,ik is the m-th column of H; Therefore at the A;-th user in the /-1h cluser, the constraints can 
be rewritten as follows: 


r H 

i.k 


hi, 


ik 


h(-i 


ik 


H-\-l,ik 


l M,ik 


o. 


Note that the dimension of H ^ is iV x (M — 1) since it is a submatrix of H, ^ formed by removing one 
column. As a result, v i)k can be obtained from the null space of H,.fc, i.e., 




( 12 ) 


’it is worth pointing out that the BS still needs to know the order of the users’ effective channel gains in order to implement NOMA as 
shown in but this imposes a much less demanding requirement compared to knowing all the users’ channel matrices at the BS. 
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where contains all the left singular vectors of H, j corresponding to zero singular values, and z tk 
is a (N — M + 1) x 1 normalized vector to be optimized later. In order to ensure the existence of v ijfc , 
N > M is assumed. 

By using the above precoding and detection matrices, the SINR for the A'-th user in the first cluster is 
given by 


SINRi K = 


H U 12„,2 

0! ltK 


121 ’ 
p 


(13) 


E/ll l V ?A- h l,lA'| 2 a?,/ + |vi,if| 

where inter-cluster interference has been removed. 

At the k-th users, 1 < k < K, the messages sij, K > j > (k + 1), will be detected with the following 
SINR: 


SINRi, = 


, r H l, | 2„,2 
v i,fcni,ifc| a iJ 


E J- 1 \^ r H t. 12 ^,,2 I I ,, 121 

/=1 l v l,fcl 1 l,lfc| a i,l + l v l,fc| p 


(14) 


If successful, Sj t i will be removed from the k- th user’s observation, and SIC will be carried out until its 
own message is decoded with the SINR, SINR\ k . 

The first user in the first cluster will decode the other users’ messages with SINR\ 1 , K > j > 2. If 
successful, it will decode its own message with the following SINR: 

I,,if h 12 2 

C r,rnl l V l,i n i,ll| “l,! 

SINR 11 = p — - - — 2 -. 


Vil 


(15) 


As can be observed from the above SINR expressions, z lk determines the SINRs through v')..!!,. 


ik 


Therefore, one possible choice of z l k can be obtained by using maximal radio combining (MRC) approach. 
Particularly, the choice of z ijk based on MRC is given by 


UKh 


^ i.k 




U"h. 


(16) 


i h ll i,ik | 


The following theorem provides an exact expression for the outage probability achieved by MIMO-NOMA 
and its high SNR approximation. 

Theorem 1. Assume that the users in each cluster are ordered as in ©. With MIMO-NOMA, the outage 
probability experienced by the k-th ordered user in the i-th cluster is given by 


, k _ ^ (-I)"*! 

Pi,t = Sv P J(K-k)\(k-l)\(K-k + p+iy 


7 (aT—M+ l,e* fc ) 
(. N-M)\ 


K-k+p -\-1 


(17) 


if Oil j > Pi, 3 > f or al1 k < j < K, otherwise P° K = 1, where e t . k = 2 Ri ’ k - 1, /3 i>k = e ijk Y^k=\ a lk> 7 (‘) 


denotes the incomplete gamma function, e* k = 


max 


£j,K 


^i.k 


p( a i, K -Pi, K ) ’ ’ p{al k -Pi,k) 


, for 2 < k < I\ and 












7 


e iX = max 


given by 


tj,K _ _ _ _ f-i, 2 6j,l 

p{ a2 i,K~Pi,K) ’ ’ P(°i,2 “A, 2) ’ P"if,l 


. A high SNR approximation for the outage probability is 


K\ 


p° ~ 


v JV-M+l ' 


(AT-M+1)! 


lf-fc+1 


( 18 ) 


(K-k)\{k- l)\(K-k + l)' 

Proof: Please refer to the appendix. ■ 

A benchmarking scheme based on conventional MIMO-OMA can be described as follows. The MIMO- 
OMA transmission consists of K time slot. During each time slot, M users, one from each cluster, are 
served simultaneously based on the same manner as described for MIMO-NOMA. As a result, the SINR 
at the A>th user in the i -th cluster is given by 




SINR; k = 


E M I H TT 


|2l 

P 


(19) 


+ Ni,k\ 

Note that the MRC detection vector used for MIMO-NOMA is also applicable to MIMO-OMA. In addition 
consider that the users in one cluster are also sorted as in ([6]). The outage probability achieved by this 
version of MIMO-OMA can be obtained in the following corollary straightforwardly by following the 
steps in the proof for Theorem Q] 

Corollary 1. Assume that the users in each cluster are ordered as in ©. By applying conventional 
MIMO-OMA, the outage probability experienced by the k-th ordered user in the i-th cluster is given by 


k -1 

p?,* = E 

p =0 


(-l)Aftl 


(. N-M)\ 


K-k-\-p -\-1 


k-r 

p )\K -k)\{k~l)\[K -k + p+1)’ 


( 20 ) 


cxKRi -1 

where (fj, = ---. A high SNR approximation for the outage probability is given by 


K\ 


P° 

r M 




,N-M+ 1' 


(n-m+i)\ 


K-k+l 


( 21 ) 


' k ~ (K — k)\(k — 1)\(K — k + 1)' 

As can be observed from Theorem Q] and Corollary [U MIMO-NOMA can achieve a diversity gain of 
(N — M + 1 )(K — k + 1)), the same as conventional MIMO-OMA. But this diversity gain is achieved 
by allowing all the K users from the same cluster to share the same bandwidth resource, which yields 
better spectral efficiency. For example, the simulation results provided in Section [VI] demonstrate that 
MIMO-NOMA can achieve a smaller outage probability compared to conventional NOMA. The superior 
spectral efficiency of MIMO-NOMA can also be demonstrated by the fact that it can realize a larger sum 
rate, as shown in the following section when the impact of user pairing is investigated. 



















IV. The Impact of User Pairing 


User pairing has the potential to reduce the complexity of NOMA systems. Specifically the users in one 
cluster can be divided into groups with fewer users in each group. A hybrid multiple access scheme can be 
used, where NOMA will be implemented among the users within each group, and conventional OMA can 
be used for inter-group multiple access. In addition to reducing system complexity, user pairing/grouping 
can also significantly increase the performance gain of NOMA over conventional MIMO-OMA, as shown 
in the following. 

In order to obtain some insightful analytical results, we focus on the case in which two users are paired 
together for performing NOMA in each cluster. Particularly the n-th and A-th ordered users from each 
cluster are scheduled to perform NOMA, where the n-th user has a better channel condition, i.e., n < k. 
By using the same choices of the precoding and detection matrices, the SNR for the k -th user in the first 
cluster is given by 


SNR ltk 

and the SNR at the n-th user is given by 


\v H 

■hi,ife 


l v U- h i,ifc| 

1 ^1 ,n 

+ 

l v l,/c 

12 1 
P 


( 22 ) 


SNR hn = p 


\ v 1 ,n 

hi, In 

l 2 a? 


|vi,„| 

2 


(23) 


conditioned on the event that the n-th user can decode the other user’s information correctly. Note that 
the power allocation coefficients satisfy a\ n + a\ k = 1. 

We are particularly interested in the sum-rate gap between MIMO-NOMA and conventional MIMO, 
which is given by 


M 

A A J2 t log ( X + SNR i,k) + log (1 + SNRi , n )] (24) 

i— 1 

i M 

- 2 [ log i 1 + p\rfk h i,ik\ 2 ) + log (1 + |vf n h ljin | 2 )] . 

i —1 

Following the same definitions used in the proof for Theorem Q] the average sum rate gap can be 
expressed as follows: 


E {A} = ME {log (1 + SNR hk ) + log (1 + SNR hn )} 

- { l0g + ^l V ffc h l,lfc| 2 ) + l0g i 1 + l V Rn h l,ln| 2 ) } 

= ME |log ^1 + ^ i j + log (l + x n a\ n p) 

M 

- y £ { log i 1 + pXk ^> + log i 1 + pXn )} ’ 





where x k = Ivf^hgifcl 2 for notational simplicity. 
After some manipulations, we can write 

M 


£ {A} = — £ {log (1 + px k )} + M£ {log (l + px n a\ n )} 
- M£ {log (l + px k a\ n )} - {log (1 + px n )} . 


( 25 ) 


The key for evaluating the rate gap £ {A} is to characterize £ {log (1 + x n (p)} which can be calculated 
as follows: 


£ {log(l + x n f)} 

/OO 

log (1 T xcj)) d(l - F Xn (x)) 
</> ri-F Xn (x) 


( 26 ) 


In 2 


l 0 1 + X(j) 


-dx. 


By applying the cumulative distribution function (CDF) of the channel gain, x n , provided in 
the proof for Theorem Q] the sum rate gap can be expressed as follows: 


£{log(l + x n cf>)} 

* f 0 X h k (*) [F ik {x)] K ~ n [1 - F ik (x)] n ~ 1 dy 


in 


(27) 


In 2 


In 2 


1 + x<f> 


1 + X(j> 


dx 


x( _ ir l3MF^F ]dx , 


where 7 n = 


K — n + p + 1 

K! —-y and the CDF F Xk (x) is obtained following the density function in 


By using 


(K—n)\(n— 1)! x k' 

the above equation and with some straightforward manipulations, the ergodic rate gap can be obtained in 
the following lemma. 

Lemma 1. Suppose that the n-th and k-th users are grouped to perform MIMO-NOMA. The average 
sum rate gap between MIMO-NOMA and conventional MIMO-OMA is given by 

M 


£ {A} = — ip(k, p) + Mp(n, pa l n ) - Mp(k, pa lt , 
M 

- y<WP), 


where 




In 2 ,/ 0 1 + x(j) 

x(-l) p - 


7(JV—Af+l,x) 
(. N-M)\ 


n— 1 / 

p=o v 1 
K—n+p+ l \ 


( 28 ) 


( 29 ) 


K—n+p+l 


dx. 
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While the analytical result in Lemma Q] can be used to replace Monte-Carlo simulations for performance 
evaluation, this is still quite complicated due to the integrals and special functions. In the following, some 
case studies will be carried out in order to obtain some insight into MIMO-NOMA. 


Case studies for the sum-rate gain of MIMO-NOMA 

In this subsection, we focus on two extreme cases as described in the following: 

• Case I: In each cluster, pair the user having the worst channel condition with the one having the best 
channel condition, i.e., n = 1 and k = K. 

• Case II: In each cluster, pair the user having the best channel condition with the one having the 
second best channel condition, i.e., n = 1 and k = 2. 

In conventional MA systems, scheduling users with better channel conditions is beneficial for improving 
system throughput, but we can show that NOMA has a behavior different from conventional MA. 


Lemma 2. For the case with N = M = 2, n = 1 and k = K, the average sum-rate gap between 
MIMO-NOMA and MIMO-OMA is given by 


S {A} = - log(e)e p E; 


fx*i i 


K 

~T~ 


2 log(e)e /> “v i E; (- \- 

V P a i,i 

N(± (fV-D'eV ‘ 


(30) 


where Ei(-) 
follows: 


denotes the exponential integral function. At high SNR, the gap can be approximated as 


£ {A} « log K + J2 logl (31) 

i=i ^ ' 

For the case with N = M = 2, n = 1 and k = 2, the average sum-rate gap between MIMO-NOMA and 
MIMO-OMA is given by (l32l) . At high SNR, the average gap can be approximated as follows: 

£{A} si K ^(-lJOogp (33) 

+ t (f)(-!)'log i). 

Proof: Please refer to the appendix. ■ 
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Define w(k) = Ym=i (j)( — l) z log/ which is a mono-increasing function of k. Lemma [2] shows that, at 
high SNR, the sum-rate gap for Case I can be approximated as (log A' + zu(K)), which means that the 
larger I\ is, the more gain MIMO-NOMA can offer compared to conventional MIMO-OMA. On the other 
hand, numerical results show that the value of k{w(k ) — w(k — 1)) quickly goes to zero by increasing k, 
which means the sum-rate gain offered by MIMO-NOMA for Case II is diminishing with increasing K. 
These two extreme cases demonstrate that careful user pairing is critical for MIMO-NOMA to outperform 
conventional MIMO-OMA. Detailed numerical analysis will be provided in Section [Vlj 

V. Cognitive Radio Inspired MIMO-NOMA 

In the previous sections, fixed choices of power allocation coefficients have been considered, and in 
this section, more sophisticated choices will be used. Without loss of generality, we focus on the same 
case as in Section HV] i.e., the n-th and A-th users from each cluster are selected to perform NOMA and 
the A-th user has poorer channel conditions, i.e., n < k. 

An important observation is that there is a dilemma in NOMA systems for choosing a 1)fc . From the 
perspective of the overall system throughput, an ideal choice of a\ t k is a^k = 0, i.e., all power is allocated 
to the user with better channel conditions. But this choice completely ignores the user fairness, and in 
this section we focus on two choices of a h k inspired by the concept of cognitive radio networks. 


A. To meet a fixed QoS requirement 

Consider that there is a targeted SINR threshold to ensure the QoS requirement at the A-th user, i.e., 
SINR i k > ej fc. This SINR requirement imposes the following constraint on the power coefficient af k \ 

i,k (|vf fc h Mfc | 2 + i) 


1 > > 


| 2 (1 + G,fc) 


(34) 
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In this paper, we will simply set a itk as follows: 


a 


i,k 


= min < 1, 


i,k (|vgfeh»,ifc| 2 + j) 


| 2 (1 T e i,fc) 


(35) 


This choice of a^ k means that the BS will give the A-th user the minimal transmission power needed to 
meet this user’s QoS requirement, and then allocate the remaining power to the n-th user. 

The outage probability experienced at the k -th user is equal to ^ = 1) or equivalently 


> 1 ) , i.e., the A'-th user’s targeted data rate cannot be supported even if the BS 


Wfj,h i4k \ 2 (l+e i:k ) 

allocates all the power to this user. Following the proof of Theorem Q] it is straightforward to show that 
a diversity order of (N — k + 1)(N — M + 1) is achievable at the k -th user, because 

i,k (|vf fe h Mfe | 2 + i) 


^i.k 


i H h , 2ri ,, , >i =p KA«r< —)■ 

The following theorem demonstrates the achievable diversity order at the n-th user. 

Lemma 3. With the cognitive radio inspired power allocation coefficient a i)k in (031) . a diversity gain of 
(N — M + 1)(K — k + 1) is achievable at the n-th user. 


Proof: Please refer to the appendix. ■ 

It is important to point out that the diversity gain at the n-th user is constrained by the k -th user’s 
channel condition due to the use of (031) . which is consistent with the finding in iflOl . Recall that cognitive 
radio inspired NOMA with signal-antenna nodes can achieve a diversity of (K — k+1) for both users IflOl . 
Therefore one advantage of MIMO-NOMA is that a larger diversity order can be achieved. In addition, 
the use of MIMO-OFDM can ensure that more users are served simultaneously. 


B. To meet a dynamic QoS constraint 


Another choice for the QoS requirement at the k -th user is to ensure the following constraint: 


log 


, , \ V l,kKlk\ 2 ^lk \ 1 , | Hu 12 \ 

1 + I ..H U 2^,2 , 1 > o l0g l 1 + l V l,fc h blfcl p) > 


(36) 


,,,JJ U 12.,2 I J 

|v life ni,ife| ol 1u + 

which means that the k -th user is willing to perform NOMA with the n-th user only if it can achieve a 
larger rate compared to the case with conventional MIMO-OMA. 

With the same notation as before, the above constraint can be expressed as follows: 

(1 + x k p) 2 


(1 + x k a\ pf 


> (1 + x k p ), 


(37) 







13 


from which the constrain for the power coefficient a\ n can be obtained as follows: 

0 < a hn < 


v/1 + PXk - 1 


px k 


(38) 


Note that 0 < ^ l+pXk 1 < 1 for arbitrary choices of x k , and therefore 

— pXk — J pXk 

we will set 


y/l+pXk-l 


®l,n 


yi + pxk - 1 
px k 


is a feasible choice. So 


(39) 


which is the maximal value of the power allocation coefficient given the constraint in (1371) . 

We first focus on the impact of this power coefficient on the outage probability at the k -th user, which 
is given by 


p ^ p|iog|i+ ^Sy |<fli ‘ 

= P (log y/lT~pXk < Rl,kj ■ 


(40) 


An important conclusion from (l40l) is that the use of the power coefficient in (l39l) ensures that the /. -th 
user experiences exactly the same outage probability as the case with conventional MIMO-OMA. This 
observation is expected since the choice of a^n is to ensure the constraint in (l36l) . i.e., the k -th user’s rate 
should not be reduced because of the use of NOMA. Following steps similar to those used in the previous 
section, it is straightforward to show that the diversity gain of this user is (N — M + 1 )(K — k + 1). 

Because the expression for the power allocation coefficient in (l39l) is very complicated, an exact 
expression for the outage probability achieved at the n-th user is difficult to find, but the achievable 
diversity gain can still be obtained as shown in the following lemma. 


Lemma 4. In the proposed CR-MIMO-NOMA system with the dynamic QoS constraint in (l36l) . a diversity 
order of (N — M + 1)(K — k + 1) is achievable by the n-th user. 


Proof: Please refer to the appendix. ■ 

It is worth pointing out that the diversity order provided in Lemma |4] is only an achievable one. After 
carrying out computer simulations, we observe that this diversity lower bound is tight for the case of 
Ri'k > 1, and a diversity order larger than (N — M + 1 ){K — k + 1) can be achieved for 0 < Ri^ < 1. 
A possible reason for this is that a loose bound is used to get the achievable diversity gain for the case 
of 0 ^ Ri,k — 1> as shown in (1871) . 
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VI. Numerical Results 

In this section computer simulation results will be used to demonstrate the performance of MIMO- 
NOMA and also verify the accuracy of the developed analytical results. For notational simplicity, we 
omit the index of the cluser, e.g., R\ k is denoted by R k . 



Fig. 1. MIMO-NOMA with a fixed set of power coefficients. M = 2, N = 3 and K — 2. a\ — \ and a§ = §. BPCU denotes bit per 
channel use. 

In Figs. Q]and|2]the performance of MIMO-NOMA with fixed power allocation coefficients is studied 
first. Particularly, Fig. Q] considers the case in which there are four users grouped into two clusters, with 
two users in each cluster. Fig. [2] considers the case in which there are three clusters, with three users 
in each cluster. All users in each cluster will participate in NOMA. Fig. Q] confirms the accuracy of 
the analytical results developed in Theorem Q] and Corollary Q] In addition this figure also demonstrates 
that MIMO-NOMA can achieve better outage performance than MIMO-OMA though both reali z e the 
same diversity gain. Fig. [2] demonstrates the accuracy of the high SNR approximation results developed 
in Theorem |T] In particular, one observation from this figure is that different users experience different 
diversity orders, which confirms the diversity order results developed in Theorem [0 

In Fig. [3] the impact of user pairing is demonstrated by using the sum-rate gap between MIMO-NOMA 
and MIMO-OMA. As can be seen from both sub-figures, the exact expression for the average sum-rate gap 
developed in Lemma Q] matches the simulation results perfectly, and the approximation result developed 
in the lemma provides a tight bound at high SNR. Comparing Fig. [3la) to Fig. |3lb), one can observe 
that the impact of K on the performance gap is much different. In Fig. [3}a), increasing the number of 
the users in each group, I \, can significantly improve the performance gap between MIMO-NOMA and 
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—□— First user, Simulation, ^=2 BPCU 
■ First user, Analytical, R^2 BPCU 
. . First user, Approximation, ^=2 BPCU 

—e—Second user, Simulation, R 2 =1.2 BPCU 
o Second user, Analytical, R 2 =1.2 BPCU 
. -o- Second user, Approximation, R 2 =1.2 BPCU 
—$—Third user, Simulation, R 3 =0.9 BPCU 
0 Third user, Analytical, R 3 =0.9 BPCU 
. . Third user, Approximation, R=0.9 BPCU 


18 20 22 
SNR in dB 


Fig. 2. MIMO-NOMA with a fixed set of power coefficients. M = 3, N = 3 and K = 3. a\ = l, = 4 and = i. 




(a) Case I: n = 1 and k = K (b) Case II: n = 1 and k = 2 

Fig. 3. The performance gap offered by MIMO-NOMA, M = 2 N = 2 and a\ = \. 


MIMO-OMA. Specifically a gain of 2 bits per channel use (BPCU) can be obtained when there are 2 
users in each group, and this gap can double when there are 5 users in each group. The reason for this 
performance gain is because we schedule the best user and the worst user, i.e., n = 1 and k = K, and the 
two selected users’ channel information becomes very different when increasing K, which is beneficial 
to the implementation of NOMA. On the other hand, Fig. Ob) demonstrates that the performance gain of 
MIMO-NOMA is diminishing with increasing K. This is because the user with the best channel conditions 
and the one with the second best channel conditions are scheduled. When increasing K, the two users’ 
channel conditions become significantly similar, which will reduce the performance gain of NOMA. 

In Fig. 0 the performance of cognitive radio (CR) inspired MIMO-NOMA for meeting the fixed QoS 
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Fig. 4. Cognitive radio inspired MIMO-NOMA with a fixed power constraint, = 1, M = 2, N = 2, K = 3 and Ri = 2 BPCU. 

requirement in (l34l) is studied. In this figure three types of curves are provided, one for P° as studied 
in Lemma [3] one for P (o 2 = 1), and one for p{N _ M+ l)( K _ k+1) . The last is provided to demonstrate the 
achievable diversity order. As can be seen in the figure, the curves for P° are parallel to the ones for 
(N-M+i)(K-k+i) > which demonstrates that the achievable diversity order obtained in Lemma [3] is tight. An 
interesting observation from the figure is that P(a 2 = 1) is a tight lower bound of P°, particularly at high 
SNR. This is because CR-MIMO-NOMA tends to satisfy the k- th user’s QoS first and therefore the event 
a 2 = 1, he., the BS allocates all the power to the fc-th user, is dominant among the three types of events 
described in the proof for Lemma [3j 

Finally, the performance of CR-MIMO-NOMA in meeting the dynamic QoS requirement in (l36l) is 
investigated in Fig. [5] Again the curves for and V are provided to facilitate the analysis of diversity 
orders. Both sub-figures demonstrate that a diversity order of (N — M + 1) ( K — k + 1) is achievable 
regardless of the choice of Rk, which confirms the accuracy of Lemma 01 Furthermore, this diversity 
order of (N — M + 1) ( I\ — k + 1) can be tight depending on the choice of R For example, in Fig. Oa), 
when Rk = 2 BPCU, the curves for the outage probability for the user with better channel conditions are 
always parallel to the ones for p ^ N _ M+ l) {K _ k+1) . In general, our carried out simulation studies reveal that 
the diversity order of (N — M + 1 ){K — k + 1) is exactly what CR-MIMO-NOMA can realize in the case 
of Rk > 1. However, in the case of 0 < Rk < 1, a diversity gain larger than (N — M + 1 )(K — k + 1) can 
be achieved, as shown in Fig. Ob). As discussed in Section [VI the reason for this is because the upper 
bound used in the proof for Lemma [4] is loose in the case of 0 < Rk < 1. 
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(a) K = 3, n = 1 and R n = R k = 2 BPCU (b) K = 2, n=l, k = K and R n =2 BPCU 

Fig. 5. CR with a dynamic power constraint, M = 2, N = 2 

VII. Conclusion 

In this paper, we have studied the application of MIMO to NOMA systems. A new design of precoding 
and detection matrices for MIMO-NOMA has been proposed, and its performance has been analyzed. 
To further improve the performance gap between MIMO-NOMA and conventional OMA, the use of 
user pairing has been considered in NOMA systems and its impact on the system performance has also 
been characterized. The cognitive radio inspired choices for power allocation coefficients have also been 
proposed to meet various QoS requirements. Simulation results have been provided to demonstrate the 
accuracy of the developed analytical results. In this paper, it is assumed that users have been randomly 
divided into multiple groups, and an important future direction is to study the design of low complexity 
approaches for dynamic clustering/grouping in MIMO-NOMA systems. 

Appendix A 
Proof for Theorem □ 

The proof can be completed in four steps. 

A. Density function of effective channel gains 

Without loss of generality, we only focus on the users in the first cluster. First recall that these users 
have been ordered according to the criterion in © which can be rewritten as follows: 


X\ > • • • > Xk 


( 41 ) 
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where Xk — |vf fc hx lfc | 2 . Define Xk as the unordered counterpart of Xk . Given the choice of v 1;fe = Ux^zx., 

and Zl - fc = iu^ fchlilfc r we have 


Uf>i, 


„ r H r, 12 _ 
v l,fc n l,lfcl — 


lfc 


= |U" hx, 


lfc | 


|U^h ui | 

An important observation is that Ux^ contains the (N — M + 1) orthogonal singular vectors, i.e., 

U^Ux, fc = IjV-M+1, 


and also note that Ux^ is independent of hu. Therefore vf^liu/. represents a unitary transformation of a 
complex Gaussian vector, which means that vf fe hx^ is still an (N — M + 1) x 1 complex Gaussian vector 
lfl3l . Therefore this unordered variable, x k , follows the chi-square distribution, and thus the probability 
density function (pdf) of x k is given by 


and its CDF is F^(x) 
pdf d: 


fx k (x) = 


(.N -M)\ 


x 


N-M 


(42) 


Jo fx k (y)dy. Therefore the ordered variable, Xk, in (HTI) follows the following 


K\fx k {x)[F ik (x)} K - k [l - F~ Xk (x)} k ~ l 
(K-k)\(k-l)\ 


(43) 


B. A unified outage probability expression 

Because the users in one cluster carry out different detection strategies, the outage probabilities achieved 
by different users will be evaluated separately first and then a unified expression for these probabilities 
will be developed. 

1) Outage probability at the user with the worst channel condition: The outage probability for the 
Ji-th user in the first cluster is given by 


P (SINR 1>K < e hK ) 


= P 


Xk&i k 


ZK22k =1 °Tfc + - p 


< t\,K 


The above outage probability can be written as follows: 


(44) 


P (SINR^k < e hK ) 


P ( XK < t * ei ’ K « > 

\ P{ a l,K-Pl,K) J 


where f3 l>K = e 1/K J2k =1 


Jf-l 2 


a 


i ,fc- 


if a l,K > Pl,K 

otherwise 


( 45 ) 
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2) Outage probability at the k-th user, 1 < k < K: The k -th user needs to decode the j -th user’s 
message, j > k, before detecting its own message. The overall outage probability for the fc-th user to 
decode its own message can be expressed as follows: 

p;, fc = 1 - P {SINR{ k > e ld , Vj e {k, ■ - ■ , K}) (46) 


= 1-P 


XkOtlj 


T > e !j5 Vj e {k, ■ ■ ■ , K} I . 


* fc ECi«?,« + P 

Following steps similar to those used in the previous subsection, the above probability can be rewritten 
as follows: 


x k alj 


x k EL I 1 a?,! + 

P ( x k > 


j > dj,Vj G {k, ■ ■ ■ , K} 

p 


(47) 


0, 


if Cl 
otherwise 




G«? 

k < j < K\. The outage probability can be expressed as follows: 
p { x k < L,k) > if C1 


where the condition, Cl, denotes of ■ > /3i j, for all k <j < K, and = dj Ejt=i a i,i 
Define e\ k = max 


po _ 
r i,fc — 


(48) 


1, otherwise 

It is interesting to observe that the expression in (l45l) is a special case of (l48l) . It is straightforward to 

evaluate that the outage probability expressions in (l48l) can also be used for the user with the best channel 

condition by letting el , = max < , L 1K a —v, ■ ■ ■ , , 9 £l,2 0 —v, -^f- 
1,1 \p{ a l,K-Pl,K) P{ a l,2-Ph2) P a il 

C. Obtaining an exact expression for the outage probability 

When the conditions, a\ k > /3i t k, are satisfied, the outage probability is given by 


fc-i 


ph =E ( k - 1 ) (- 1 )” 

p =0 

f e l,k 


K\ 


x 


Jo 
k -1 


p j ■ (K - k)\(k - 1)! 
fx k ( x ) \ F x k ( x)} K ~ k+p dx 


(49) 


^,k- 1\ (-1 yK\[F, k (e\ jk )] 

p j(K-k)\(k-i)\(K-k+p+iy 


P 

By applying the CDF of the unsorted variable x k , we obtain 

_ “ ( k - A (~ 1 ) p k\ hMdv 


p =0 


K-k-\-p -\-1 


p J (K - k)\(k - 1)\(K - k + p + 1) 


By applying the incomplete gamma function, the exact expression of the outage probability can be obtained 
as in the theorem. 
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P° 

r l./ 


k -1 

E 

p =o 

fe-i 

E 

p =o 


k - 1 
P 


k- 1 
P 


-1)AA! 


(JV - M)\ ( 1 - e-M 


- X—fc+P+1 


(AT - fc)!(fc — l)\(K — k + p + 1)((N — M)\) K ~ k+ P +1 


(50) 


-1) P A3 — 




MX) 


.0' 


X— k-\-p-\-l 


(.A - A:)!( fc - 1)K A ' - k + p + 1)((N - M)\)) K ~ k +P+ 1 


D. High SNR approximations 

By applying the series expansion of the incomplete gamma function llT5l . the outage probability can 
be first expanded as in (l50l) . At high SNR, the outage probability can be approximated as follows: 


fc-i 


p?,*=E (\ K- 1 )^ 1 

p=0 ' ^ 


x 


p-q.fe v°° 
e l^q=N- 




q=N-M+l 


K—k+p+l 


(K - k)\(k - 1)\(K - k + p+ 1)((N - M)\) K ~ k+ P +l 

Ef; 1 ) (- 1 )'" 

p= o \ ^ / 


X 


(AT-M)!( e j |fc ) 


N-M+l 


(JV-M+1)! 


K—k+p+l 


(.A - k)\(k - 1)!(AT -k + p+l)((N- M)\) K ~ k +P + 1 


A'! 




JV-M+1 


(JV-M+1)! 


A'-fc+l 


(A" — k)\(k — 1)\(K — k + 1) 
Therefore the theorem is proved. 


(51) 


Appendix B 
Proof for Lemma [2] 

The sum-rate gap will be evaluated separately for two cases in the following subsections. 


A. Case I with n — 1 and k = I\ 

First recall the following integral from Eq. (3.352.4) in 1(1511 : 



1 

1 + X(j) 


e lx dx 




l 

~4> 
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By using the above result and also Lemma [H <^>(1,0) can be expressed as follows: 




In 2 ,/ 0 1 + xcj) 

K 

E 


1 — 7i 


[7(1,®)] 


K 


In 2 


. 1=1 




K 


l 


dx 


(52) 


For the worst user, a direct use of Lemma [Q results in a quite complicated expression for ip(K,<f>). 
Instead, we can find a simpler alternative way to calculate this factor, as shown in the following: 


<p(K, <t>) = / log (1 + x<f>) f XK {x)dx 

Jo 

poo 

—K log(e) / In (1 + x(j>) e~ Kx dx 

Jo 


(53) 


= - log(e)e * Ei - 


K 


By substituting (l52t and (l53l) into the expression for the rate gap, the expression in (l30l) can be obtained. 
To obtain the high SNR approximation, first recall that the exponential integral function has the following 
series representation 03: 


Ei(x) = C + ln(-x) + V 


3 = 1 


for x < 0. Therefore at high SNR, we have the following approximation: 

( K' 


£ {A} ~ — log(e) ( C + In ( 

K 


V P 


(54) 


2 

M2 


E 


. 1=1 

2 log(e) ( C + In 


-1)' C + In 


l 


Puli 


K 


1^2 (E 


1=1 


K 


I< 

pu\i 


-l)‘(C + ln(i 
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After some manipulations the average gap is given by 

^+~C + ln(+)-2>„< 

log e V P J 

+ E(TM‘fc + in" 


1=1 

K 


1=1 


K 


E (A (-i )' hii + 2 E (A (-i)'(- to “h) 


i=i 


K 


= C + ln| y I - 21na?, + + ( K ll-li'ln/ 


'u + E (i) (- 1 ) 

i=i v 7 


C + In ( - 

P 


K 

E 


1=0 


K 


-l)‘(C + ln(- 


K 


1=0 


K 


2(-i„ a y + 2E 


After removing some common factors, the average gap can be simplified as follows: 

£{A} - C + In f-V 2111a?,- (c + h/ 1 


loge 


V P 


P 


K 


1=1 


K 


( / (-l) / ln/ +2111a?! 


K 


1=1 


K 


InA' + E , (- 1 )' 1 "'- 


And the first part of the lemma is proved. 


(55) 


(56) 


B. Case II with n = 1 and k = 2 

It is more complicated to evaluate the average gap for Case II due to the complicated expression for 
x 2 • In particular, the factor ip(2,cj)) can be expressed as follows: 

POO ^ 


+(2,</>) = 

+72 


In 2 J 0 1 + xcj) 
[1 - e~ x ] K 


1-72 


[1 - e~ x ] K ~ l 

K - 1 


K 

POO ^ 


dx 


K-l 


In 2 J 0 1 + x(j> 


i -*E 


p=0 


Ji -1 


p 


-l) p e 


—px 




(57) 
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Since JT TT^ dx —> oo, therefore it is important to remove the factors in the integral related to y-h- in 
order to facilitate the high SNR approximation. Motivated by this, the factor <p(2, (ft) can be rewritten as 


follows: 


ip{2,cft) 



(58) 


Following the steps as those used in the previous section, the integral in the above equation can be 
evaluated and we can have the following: 


ip(2,(ft) 



(59) 


Substituting (l59l) and (l52l) into (|28T) . the exact exact expression of the average rate gap can be obtained 
as in the lemma. 

At high SNR, the exponential integral function can be simplified as discussed previously, and the 
average rate gap can be approximated as follows: 



With some algebraic manipulations, the average rate gap can be expressed as (I6TI) . After those common 
factors in (I6TT) are removed, the high SNR approximation shown in (l33l) can be obtained, and the lemma 
is proved. 
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Appendix C 
Proof for Lemma [3] 


Without loss of generality, we will take the users in the first cluster as an example. Recall that the n-th 
user, n < k, in the first cluster tries to decode the A:-th user’s message with the following SINR: 


SINR k hn = 


r H u 12^2 
grp 1 !,Ira | a l,k 


H F 12^,2 , I 

l,n n l,ln| a l,n + ~ p 


1 • 

P 


If successful, the n-th user will decode its own message with the following SINR: 


SINRl n = p|v^ n H lilPl | 2 ay 


12 „ 2 


= p| V l,n H l,lPll 


1 - 


ei ,k (|vff fc hi ) i fc | 2 + i)\ 
|vff fc h lilfc | 2 (l+ e lifc ) J' 


(62) 


if a 2 n > 0. Again with the same notation used in the proof of Theorem [H the SINR at the n-th user can 
be expressed as follows: 


SINRl n = px n 




1 - 


X k + 


Xk( 1 + ei ifc ) 


(63) 


if a\ n > 0. 

Thus there are two conditions before the SINR expression in (1631) can be used. One is o 2 n > 0 and the 

other is that the n-th user can decode the k -th user’s message, i.e., log(l + SINR \ ,n) > /21 /, . Therefore 

the outage events at the n-th user can be categorized into three following types: 

1) Events with a\ n = 0, which means 1 < or equivalently 

^ e L 
Xk < — 

P 


2) Events with a\ n > 0 and log(l + SINR^ n ) < Ri )k . 

3) Events with a\ n > 0, log(l + SINR^ n ) > and SINR^ n < ei t7l . 


( 64 ) 
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Because |vf fc hi ilfc | 2 < |v^ n h ljln | 2 , it is straightforward to show SINR\ n > SINR\ k , which means 

P (al n > 0, log(l + SINRlJ < R 1>k ) (65) 

= E 0<aln <i{P (SINR* tn <e ltk )}) 

= £ 0<a i n <i {P (SINR* >k = e 1>k , SINR k hn < e 1>fe )} = 0, 

i.e., the n-th user can decode the k -th user’s information as long as the k -th user can decode its own. But 
if a\ k = 1, i.e., the BS allocates all the power to the A-th user, outage will occur at the r/-th user. 
Therefore the outage probability experienced by the n-th user is given by 

€l,k (xk + ^ 


K = ? PXr, 


Xk( 1 + ei.fc) 


. _ £ i ,k 

^ ^l,n? %k ^ 

P 


p ( X k < — 

\ p 


( 66 ) 


which follows from the following simplification: 

P (a 2 hn > 0. log(l + SINRlJ > < (,„) 

= P K„ > 0, SINRl„ < £l ,„) . 

Define the first factor in the expression for the outage probability in 
) < e i ,n,Xk > -y- I • This factor can be evaluated as follows: 


(67) 


by Q 2 = 


P PX r 


Xk 


Q'2 — P ( X r 


6l k 

Xk 


^ ^l.nj %k ^ 


£l ,k 

P 


^Xk i P I %k ^ %n < ^ 


£l,n%k 


where i ln — 


range 


(l+£l,fc)£l.n 


. It is important to note that the expectation in 


£ i,k . . £ i.k . ~ 

- < X k < -h £i n , 

P P 


( 68 ) 

(69) 

is taken over the following 


where the upper bound is due to the constraint x k < 
Now the factor Q 2 can be evaluated as follows: 

£ l,n,X k 


£1 ,nXk 


a*-- 


£ 1 ,k 


Q 2 — £x k \ R 


X k 


e l,fc 

p 


F Xn (x k ) 


(70) 


"+Cl,rj 


£ l,nV 


£ l ,fc 
p 


y - 


£l ,k 


F xjy) f Xk {y)dy. 
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At high SNR, Q 2 can be upper bounded as follow: 


-+ei • 


since F. r 


Cl ,nV 


y-- 


mr)-F Xn (y)<F x 


Q 2 < 


gl ,nV 

l l ,k 

v ir 


e l,fc 

P 


fx k (y)dy, 


< 1. 


With this upper bound, the overall outage probability can be upper bounded as follows: 

Cl,A: 


K = Q 2 + P [x k < 


P 


< 


e l ,fc 


f e l,fc 
P 


= F, 


+ 61,73 


Cl,fc 

P 


/+?/)%+ P ( Xfc < — 


+ Cl- 


To find a high SNR approximation of F Xn {x), first recall that the CDF of x n is given by 

n—l / -■ \ rn / \-<X—n+j+l 


Fx n {%) — 7 n y ^ 

j =0 


n — l 

7 


■i) 




K — n + j + 1 


When x —y 0, we have 


4W = 


7(1V — M + 1, x) 


x 


N—M+l 


(n — m)\ (N-M + iy: 

Therefore, the CDF of x n can be approximated as follows: 


n—l 


F x n (x) ~ 

j =0 

In 


n — l 

j 


-1 y 


x n - m +1 
(7V-M+1)! 


K-k+j -\-1 




FT — n + j + 1 

JV-M+l \ K-k+1 


(71) 


(72) 


(73) 


(74) 


(75) 


K-n + 1 \(N — M + 1)! / 

when x —y 0. 

Substituting the above approximation into (1721) . the overall outage probability can be upper bounded as 
follows: 


/ 


P° < 


7n 


->■ 


K — 77- T 1 
1 


e l,fc 1 2 

T - + e D 


N—M+l \ -fC-fc+1 


(iV-Af+1)! 


(76) 


(N-M+l)(K-k+l) ' 


P 


And the proof is completed. 
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Appendix D 
Proof for Lemma [4] 

With the same notations used in the proof for Theorem [0 the SINR at the n-th user can be expressed 
as follows: 

\/l + px k - 1 


SINRl n = px r 

if it can decode the &;-th user’s message, i.e., 


px k 


(77) 


, Xn&i k 

!°g 1 + — 2 —rr > R pk- 

Xn^ln + - ' 

Therefore the outage events at the n-th user can be categorized into the two following types: 

• Events in which the n-th user cannot decode the k -th user, i.e., 

log (l+ <Ri, k 

\ X n(-4.n +pj 

• Events in which the n-th user can decode the k -th user, but cannot decode its own, i.e., 

log ( 1 + Xn ^’ k , ] > R ltk , 


(78) 


x n<Xl, n +j 


and 


log (l + px n a 2 ln < Ri tU ) . 

Therefore the outage probability experienced by the n-th user is given by 


K = P[ log 


+ P log 


X n P + 1 

^n“l ,nP+l / 

x n p+ 1 

x n a 2 ln p+ 1 


< R 


l,k 


( 79 ) 


> R ltk , log (l + px n a\ n < . 


Q 3 


The first factor of P° can be calculated as follows: 


Qa = P ( log ( Xn P + 1 ^ j < Ri k 


X nalnP +1 


T-» ( <2 R l,k— 1 

P [ Xk < x n < -—? If—- 

1, 


(80) 


if 1 > 2 Rl ’ k ai :n 

otherwise 


The constraint of 1 > 2 Rl • k ai >n is equivalent to the following one: 

- 2 ) 


x k > 


P 


(81) 
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As a result, P (log 


x X a{ + p +1 ) < k ) can be expressed as follows: 


P ( 


Xf~ ^ Xji 

Xk > 


2 i.fc —1 


p(l-2" 1 » fc ai |n ) 

2 «i,fc(2 a i,fc-2A 


= < 


( 0 W l,fe _ 1 

P Zfc < X n < , , ~~7% 


7T7T — , if 1 > 

& -Ri,fe > 1 
) , if 1 > 2 Rl ’ k ai >r 
& Rl,k < 1 

otherwise 


p(l-2" 1 .'=ai i „) 


( 82 ) 


Consequently P (log ( x ^ +1 


< 


Ri,kj can be upper bounded as follows: 


/ 2 Rl - k - 1 
Qa ^ P St < „ d t 

V p(l - 2 Hl - k a ltn ) 

2 R hk(2 R i* - 2)\ 


(83) 


+ P [x k < 


P 


if R\ t k > 1, otherwise 


/ 2 Rl ’ k - 1 

Q 4 = P \ Xu < x n < — - 5 -- 

1 p(l-2 R ^a hn ) 


(84) 


to 


( 2 Rl ' k (2 Rl > k 2) 

Comparing (1831) to (1841) . one can observe that the probability P ( x k < — ^ j does not need 
be taken into consideration for the case of 0 < < 1. In the following, we first focus on the case 

Rl,k > 1- 

It is important to note that the constraint x k < x n < ——— yields the following additional 
constraint for 


x k < 


2 Rl ’ k - 1 


which leads to the following inequality: 


oil — 2 Rl nASeLiIA 

1 V pXk J 


(85) 


x k < 


2 2R 1 ,k _ X 

p 


( 86 ) 


Therefore P (log ( x ^ P+ p+i ) < can be upper bounded as follows: 


P (log 


X n P+ 1 
Xn^lnP +1 


< R 


1 ,k 


( 87 ) 


2 R i,k - i 


< P \X k < X n < — --jj--- 

p(l - 2"i. fc «i jn ) 

2 R i,k(2 R hk - 2) 


,x k < 


2 2R 1 ,k _ l 

p 


+ P [x k < 


< P ( Xk < 


p 

2 2R l,k _ i 

P 


+ P ( Xk < 


2 R i,k(2 R hk - 2)' 


P 
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Following steps similar to those used in the previous section, we have the following asymptotic result: 

1 


2 2i h,fc — X' 

P [ X k < --- I -A 


p 


(. N-M+l)(K-k+l ) • 


( 88 ) 


The other probabilities have the same asymptotic behavior, and therefore by combining (1871) and (1881) . we 
have 

X n P + 1 \ „ \ 1 


P log 


x n al n p+ 1 


< Rl,k < 


p (N-M+l)(K-k+l ) ’ 


( 89 ) 


where a<b denotes lim > lim llT6ll . The above conclusion is also valid for the case 0 < Ri k < 1. 

p—too 10 SP p—too 10 ®/ 3 


The factor Q 3 can be upper bounded as follows: 


Qz < P (log (l + px n ot 2 ln < Ri, n )) 

2 Rl ’ n - 1' 


= P [ x k < x n < 


pain 


A hidden constraint on x k due to x k < x n < 2 ^ 1 is 




x k < 


2 r '-' 


P 


s/l+px k -l 

pXk 


which yields x k < 


2 2 ^l,n 1 


Therefore the factor Q 3 can be further upper bounded as follows: 

/ 2 Rl ’ n — 1 2 2Rl ’ n — 1 

<?3 < ? ( X k < X n < -j-, X k < - 

P 


P a l,n 

,n _ 

< P ( X k < --- | -A 


,(N—M+l)(K—k+l) 


(90) 


(91) 


(92) 


Combing 


P' 

and (l93l) . the overall outage probability can be upper bounded as follows: 
P° 


(N-M+l)(K-k+l) 


P' 


(93) 


And the proof is completed. 
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